Circular dielectric cavity and its deformations 
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The construction of perturbation series for slightly deformed dielectric circular cavity is discussed 
in details. The obtained formulae are checked on the example of cut disks. A good agreement is 
found with direct numerical simulations and far-field experiments. 
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I. INTRODUCTION 

Dielectric micro-cavities are now widely used as micro- 
resonators and micro-lasers in different physical, chemi- 
cal and biological applications (see e.g. [|J and ref- 
erences therein). The principal object of these studies 
is the optical emission from thin dielectric micro-cavities 
of different shapes [||. Schematically such cavity can be 
represented as a cylinder whose height is small in com- 
parison with its transverse dimensions (sec Fig. [T]) . If the 
refractive index of the cavity is n\ and the cavity is sur- 
rounded by a material with the refractive index n 2 < n\ 
(we assume that the permeabilities in both media are the 
same) the time-independent Maxwell's equations take the 
form (see e.g. 



V ■ B } = 







V x Bj ■ — -irijkEj , 



Vx£; = ikB, 



(1) 



where the subscript j = 1 (resp. j = 2) denotes points 
inside (resp. outside) the cavity and k is the wave vector 
in the vacuum. These equations have to be completed 
by the boundary conditions which follow from the conti- 
nuity of normal B v and nE v and tangential E T and B T 




FIG. 1: Schematic representation of a dielectric cavity. 
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components 

n\E ly — n\E 2v , B\ v = B 2v , E\ T = E 2t , B\ T = B 2t ■ 

In the true cylindrical geometry, the z-dependence of 
electromagnetic fields is pure exponential: ~ e lqz . Then 
the above Maxwell equations can be reduced to the 2- 
dimensional Helmholtz equations for the electric field, 
Ej z , and the magnetic field Bj Z along the axe of the 
cylinder 

(A + n 2 fc 2 ) E jz (x,y) = , (A + n 2 fc 2 ) B jz (x,y) = (2) 
with the following boundary conditions 

8E lz 8E 2z dB lz dB 2z 



E\ z — E 2z ,Bi z — B 



22. 



dr 



8t ' dr 



dr 



and 



1 dB lz 1 8B 22 



dE z 



dv 



dv 



n\ dE lz n\ dE 2z 



^2 



dv 



dv 



kn\n\ dr 
q{n\ - n\) dB z 
kn\h\ dr 



(3) 



Here n 2 = ?i 2 — q 2 /k 2 plays the role of the effective two- 
dimensional (in the x — y plane) refractive index. 

When fields are independent on z (i.e. q = 0) bound- 
ary conditions © do not mix B z and E z and the two 
polarizations are decoupled. They are called transverse 
electric (TE) field when E z — and transverse magnetic 
(TM) field when B z = 0. Both cases are described by 
the scalar equations 



(A + njfc 2 ) 9j(x,y) = 



(4) 



where ^(x,y) stands for electric (TM) or magnetic (TE) 
fields with the following conditions on the interface be- 
tween both media: ^1 = ^2 and 



dv 



1 a*i 

n 2 dv 



dv 

1 9*2 

n 2 , dv 



for TM polarization , (5) 
for TE polarization . (6) 
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These equations are, strictly speaking, valid only for an 
infinite cylinder but they are widely used for a thin dielec- 
tric cavities by introducing the effective refractive index 
corresponding to the propagation of confined modes in 
the bulk of the cavity (see e.g. @). In practice, it re- 
duces to small changes in the refractive indices (which 
nevertheless is of importance for careful comparison with 
experiment Q). For simplicity we will consider below 
two-dimensional equations (HJ as the exact ones. 

Only in very limited cases, these equations can be 
solved analytically. The most known case is the circular 
cavity (the disk) where variables are separated in polar 
coordinates. For other cavity shapes tedious numerical 
simulations are necessary. 

The purpose of this paper is to develop perturbation 
series for quasi-stationary spectrum and corresponding 
wave functions for general cavities which are small de- 
formations of the disk. The obtained formulae are valid 
when an expansion parameter is small enough. The sim- 
plicity, the generality, and the physical transparency of 
the results make such approach of importance for tech- 
nological and experimental applications. 

The plan of the paper is the following. In SectionHTlthc 
calculation of quasi-stationary states for a circular cavity 
is reviewed for completeness. Special attention is given 
to certain properties rarely mentioned in the literature. 
The construction of perturbation series for eigenvalues 
and eigenfunctions of small perturbations of circular cav- 
ity boundary is discussed in Section lllll The conditions 
of applicability of perturbation expansions are discussed 
in Section HVl The obtained general formulae are then 
applied to the case of cut disks in Section [V] Some tech- 
nical details are collected in Appendices. 



II. DIELECTRIC DISK 

Let us consider a two dimensional circular cavity of 
radius R made of a material with n > 1 refractive in- 
dex. The region outside the cavity is assumed to be the 
air with a refractive index of one. The two-dimensional 
equations ^ for this cavity are 



(A + n 2 fc 2 )* 
(A + fc 2 )* 



when r < R , 

when r > R . (7) 



There is no true bound states for dielectric cavities. The 
physical origin of the existence of long lived quasi-bound 
states is the total internal reflection of rays with the in- 
cidence angle bigger than the critical angle 



. 1 

arcsm — . 
n 



(8) 



To investigate quasi-bound states one imposes outgoing 
boundary condition at infinity, namely, we require that 
far from the cavity there exist only outgoing waves 



In cylindrical coordinates (r, 6), the general form of the 
solutions is the following 

iT/ , / a m J m (nkr)e™ e , r<R, 

where m = 0, 1, . . . is an integer (the azimuthal quan- 
tum number) related to the orbital momentum. J m (x) 
(vesp.Hm\x)) stands for the Bessel function (resp. the 
Hankel function of the first kind) of order m. Due to 
rotational symmetry, eigenvalues with ra/0 are doubly 
degenerated. 

By imposing the boundary conditions ([5]) or ([6]) one 
gets the quantization condition 



V J rn 



(nkR) 



(i)' 



H, 



where 



1 for TM polarization 
n 2 for TE polarization 



(10) 



(11) 



The quasi-stationary eigenvalues of this problem depend 
on azimuthal quantum number m and on other quantum 
number p related with radial momentum: k = k mp . They 
are complex numbers 



k 



iki 



(12) 



where k r determines the position of a resonance and ki < 
is related with its lifetime. 

In Fig. [2] we plot solutions of Eq. (fT0|) obtained numeri- 
cally for a cylindrical cavity with refractive index n = 1.5. 
Points are organized in families corresponding to differ- 
ent values of radial quantum number p. The dotted line 
in these figures indicates the classical lifetime of modes 
with fixed m and k — * oo. Physically these modes cor- 
respond to waves propagating along the diameter whose 
lifetime is given by 



Im(kR) = — In 



n-l 
n+1 



-0.53648 



(13) 



In the semiclassical limit and for Im(fci?) <C Re(fci?) sim- 
ple approximate formulae can be obtained from standard 
approximation of Bessel and Hankel functions [7J : 

• when m < z 



(14) 



ty(x ) tx e ife|S| when 



oo 




(15) 
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a) b) 

FIG. 2: Quasi-stationary eigenvalues for a circular cavity with n — 1.5. a) TM polarization, b) TE polarization. Filled circles 
are deduced from direct numerical resolution of Eq. 1|10[1 . while open squares indicate semiclassical approximation for these 
eigenvalues based on Eqs. (I16p and (|17p when \kR\ < m < n\kR\. 



Denoting x = Re(fci?) and y = Im(A:i?) and assuming 
that !/ C i, i » 1, and m/n < x < m one gets (see e.g. 
Q) that the real part of (fT0|) can be transformed to the 
following form 



\j n 2 x 2 — m 2 



arctan v 



m arccos ■ 



7T 
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+ (p-1)tt (16) 



where the integer p — 1, 2, ... is the radial quantum num- 
ber and v is defined in pip . The imaginary part of equa- 
tion ([TOl) is then reduced to 



(17) 



vx{n 2 -l)\H£\x)\ 



where £ = 1 for TM waves and £ = n x /(m (n + 
1) — n 2 x 2 ) for TE waves. When x and m are large, y is 



exponentially small as it follows from (JTf 

The above equations can not be applied for the most 
confined levels (similar to the "whispering gallery" modes 
in closed billiards) for which nx is close to m. In 
Appendix it is shown that real part of such quasi- 
stationary eigenvalues with 0(m~ 1 ) precision is given by 
the following expression 



rn 
n 



n 



/TON 

I 2" ) 



1/3 



v\J n 2 



20n \m) 



n 2 r\ v 



2v{n 2 - l) 3 / ; 



■( — 

\3v 2 



hi* 



(18) 

2/:S 




40 60 

Re(kR) 



100 



FIG. 3: Additional branches of quasi-stationary eigenvalues 
for a circular cavity with n = 1.5. Black circles indicate 
the TM modes and red squares show the position of the TE 
modes. Solid black and dashed red lines represent the asymp- 
totic result (|19p for respectively the TM and the TE modes. 



where r\ v is the modulus of the p^ 1 zero of the Airy func- 
tion (K2\i . 

A more careful study of (fTU|) reveals that there ex- 
ist other branches of eigenvalues with large imaginary 
part not visible in Fig. [5J Some of them are indicated 
in Fig. [3J These states can be called external whispering 
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FIG. 4: Example of a deformed circular cavity. Shaded areas 
represent regions where the refractive index differs from the 
one of the circular cavity. 



gallery modes as their wave functions are practically zero 
inside the circle. So they are of minor importance for our 
purposes. They can also be identified with above-barrier 
resonances. In Appendix [A] it is shown that in semiclas- 
sical limit these states are related with complex zeros of 
the Hankcl functions and they are well described asymp- 
totically (with C(m _1 ) error) as follows 



u m,p 



?y / \e-*-/» 



36-4-/3^2 ( 2 ^ 1/3 



\VT] V 



ip 

20 \ ;// 

c -2i7r/3 



2 ( n 2 _ !)3/a 



l-^ 2 ) (- 

3 I \m 



2/3 



(19) 



with the same r\ v as in (|18p . 

Similar equations have been obtained in Q . 



will be demonstrated below) that the criterion of appli- 
cability of the quantum perturbation theory is, roughly, 



5a k 2 < 1 



(22) 



where 5a is the area where perturbation 'potential' 5n 2 
is non zero (represented by dashed regions in Fig.0]). 

To construct the perturbation series for the quasi- 
stationary states, we use two complementary methods. 
In Section UlI Al we adapt the method proposed in fllllriT] 
for diffraction problems. The main idea of this method 
is to impose the required boundary conditions (|5|) or ([6]) 
not along the true boundary of the cavity but on the cir- 
cle r — R. Under the assumption (|22[) this task can be 
achieved by perturbation series in A. In Section [ill Bl we 
use a more standard method based on the direct pertur- 
bation solution of the required equations using the Green 
function of the circular dielectric cavity. Both methods 
lead to the same series but they stress different points 
and may be useful in different situations. 

For clarity we consider only the TM polarization where 
the field and its normal derivative are continuous on the 
dielectric interface. For the TE polarization the calcu- 
lations are more tedious but follow the same steps. To 
simplify the discussion we assume that the deformation 
function f(9) is symmetric: f(—6) = f(0) (as in Fig. 0|. 
In this case the quasi-stationary eigenfunctions are cither 
symmetric or antisymmetric with respect to this inver- 
sion. Then in polar coordinates, they can be expanded 
either in cos(p#) or sin(p0) series. The general case of 
non-symmetric cavities is analogous to the case of degen- 
erate perturbation series and can be treated correspond- 
ingly. 



A. Boundary shift 



III. PERTURBATION TREATMENT OF 
DEFORMED CIRCULAR CAVITIES 



The condition of continuity of the wave function at the 
dielectric interface states 



In the previous Section we have considered the case of a 
dielectric circular cavity. It is one of the rare cases of inte- 
grable dielectric cavities in two dimensions. The purpose 
of this Section is to develop a perturbation treatment for 
a general cavity shape which is a small deformation of the 
circle (see Fig|3}. We consider a cavity which boundary 
is defined as 



r = R + Xf(9) 



(20) 



in the polar coordinates {r,0). Here A is a formal small 
parameter aiming at arranging perturbation series. 

Our main assumption is that the deformation function 
Xf(9) is small 



|A/(0)|«J2 



(21) 



Of course, for the quantum mechanical perturbation the- 
ory this condition is not enough. It is quite natural (and 



9 1 (R + \f(6),6) = 9 2 (R + \f(6),0) 



(23) 



where subscripts 1 and 2 refer respectively to wave func- 
tion inside and outside the cavity. Expanding formally 
*i,2 into powers of A one gets 



[*i 
1 



§A 2 / 2 W 



V 2 ]{R,9) = -\f(0) 
5 2 *i 5 2 * 



dr 2 



dr 2 



5*i 5*2 
dr 

(R,0 



dr 



(24) 



For the TM polarization the conditions ([5]) imply that the 
derivatives of the wave functions inside and outside the 
cavity along any direction are the same. Choosing the 
radial direction, one gets the second boundary condition 



dr 



(R + Xf(9),9) 



5* 5 
dr 



(R + \f(0),6) 



(25) 
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which can be expanded over A as follows 

"0*1 0* 2 " 



dr dr 



(iZ,i 



A/(0) 



dr 2 



di 



(26) 



a 3 *i 5 3 * 2 

0r 3 0r 3 



We find it convenient to look for the solutions of Eqs. ([24 
and (|26[) in the following form 



*i(r, 



J m (nkr) 
J m {nx) 

J p (nkr) 



a p 



p^rn 



cos(m9) 

cos(p8) 



J p (nx) 



(27) 



008(7716*) 



M^T^cos(^). (28) 



Here and for all which follows, x stands for kR. These 
expressions correspond to symmetric eigenf unctions. For 
antisymmetric functions all cos(. . .) have to be substi- 
tuted by sin(. . .). 

From (|24| and (|26| one concludes that the unknown 
coefficients a p , and b p have the following expansions 



a p — Xa p 



X 2 (3 P 



b p = A 2 7 P 



(29) 



Correspondingly, the quasi-stationary eigenvalue, kR 
x, can be represented as the following series 



x = Xq + Xxi + X 2 X2 + 



(30) 



Here xo is the complex solution of (|10|) which we rewrite 
in the form 



S m (x )=0 (31) 
introducing for a further use the notation for all m and x 



JL 



H, 



(i)' 



S m (x) = n~p-{nx) 77r( x ) ■ 

Jrr, ffW 



(32) 



The explicit construction of these perturbation series is 
presented in Appendix [Bl The results are the following. 
The perturbed eigenvalue (f30|) is 



X = .T 



1 — XA mm + A 



y 2 O^mm B mm ) 



H 



(1)' 



+ x o(A mm — B mm ) — ^jy(xo) 



(33) 



- (n 2 - l)x ^ A r,ik 



Sk(xo) 



A 



km 



0{X 3 



The coefficients of quasi-stationary eigenfunction (|2~7| 
and J28l) are 



Xxq(ti 



S p (xq) 



~\~ A ^ Ap rn A 7n rn ^ 



1 / f Hg>' 
-B pm \l + x Q \-—j 



H, 



H, 



(i) 



Xq dS p 

S p dx 



+ xo(n 2 - 1) 2J -4p/c 



k^m 



Sk(x ) 



C(A 3 



(34) 



and 



the Wronskian relation (see e.g. Q 7.11.34) 



bp = X 2 -xl(n 2 - l)B pm 



+ C(A 3 



(35) 



H$\x)hV\x)-HW{x)HV>(x) 



4i 



(36) 



In these formulae A mn and B mn are the Fourier harmon- 
ics of the perturbation function f{9) and its square, given 
by (|B5|1 and (|B10p respectively. The above expressions 
are quite similar to usual perturbation series and S p (xq) 
plays the role of the energy denominator. 

It is instructive to calculate the imaginary part of the 
perturbed level from the knowledge of the first order 
terms only. Assuming that Im (xo) <C Re (xo) and using 



one gets from (|3"3"j) 



Im(x) = Im(a;o)(l + p) - X' 



2x (n 2 - 1) 



E 



A 2 



Ml 



where (assuming that m > xq so Hm / Hm^ is real cf 



(37) 
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m) 



— A.A mm + A 
x o{A 2 mm - B 



jf^m - B mm ) (38) 



1 ' tm ) (IT ) 

Bra 



Expression (|37[) without the p correction (which is multi- 
plied by Im(xo) and so negligible for well-confined levels) 
can be independently calculated from the following gen- 
eral considerations. From ([7]) it follows that inside the 
cavity 



n 2 (k 2 -k* 2 ) / |*(a;)| 2 df=J 



where 



J 



d d 
**tt^* - 4" — ** )da 
ov ov 



(39) 



(40) 



In the first integral the integration is performed over the 
volume of the cavity, V, while the second integral is taken 
over the boundary of the cavity, B. v is the coordinate 
normal to the cavity boundary and J represents the cur- 
rent through the boundary. For the TM polarization, it 
equals the current at infinity. 

From (HU) and ([34]) it follows that the field outside the 
cavity in the first order of the perturbative expansion is 



*2(r, 0) = — -y — -cos(m6>) 



H%\x) 



(41) 



X(n 2 



1)X A vm Hf ' y COs(pg) 

S p {x)B p '(x) 



p^m 



The current can be directly calculated from (|36p or from 
the asymptotic of the Hankel functions. Then the current 
can be written 



J = 4i 



1 



(i) 



pm 



(42) 



To calculate the integral over the cavity volume in the 
leading order, the non-perturbed function can be used 
inside the cavity. Then the integration over the circle 
r = R leads to 



v J™(nx ) J 



J^(nkr)rdr . (43) 



The last integral is (see [7j 7.14.1) 
i-R 

J^(nfcr)rdr 



R 2 



J£(nkR) + J^(nkR)(l 



{nkRf 



From the eigenvalue equation (flQ|) and the asymptotic 
([To]) , it follows that 



' J m (nkR) 



Therefore 



J 2 n (nkr)rdr « J^(nxo 



(n 2 - l)i? 2 



o 



2n 2 



Combining these equations leads to 



Im x = 

7T 



1 



(n 2 - ljarol^, 



(i). 



(44) 



(45) 



AV-l)z £ 



-4 2 



Snff'i 



Wi 



According to (jTTJ) the first term of this expression is 
the imaginary part of the unperturbed quasi-stationary 
eigenvalue (assuming that Im (xq) <C Re (xo)) and one 
gets (|37[) using only the first order corrections. The miss- 
ing terms are proportional to the imaginary part of the 
unperturbed level and can safely be neglected for the well 
confined levels. 

These calculations clearly demonstrate that the defor- 
mation of the cavity leads to the scattering of the ini- 
tial well confined wave function with Re(/ci?) < m < 
nKe(kR) into all possible states with different p mo- 
menta. Among these states, some are very little confined 
or not confined at all. These states with p <Re(kR) give 
the dominant contribution to the lifetime of perturbation 
eigenstates. Such scattering picture becomes more clear 
in the Green function approach discussed in the Section 

mm 

The important quantity for applications is the far-field 
emission. It is calculated using the coefficients a p (|3"4")1 
and bp ((33)) in (f2"8"|) and substituting its asymptotic for 

H^ikr): 



Then one gets 
where 



irkr 



irkr 



3 i(fer— -rrp/2— -it/4) 



j(kr-ir/4) 



F{6) 



(46) 



(47) 



p— YKmjl 

F{9) = (i + b m )_ cos(m0) 

+ LK + M (1) COs(pg) 
p^m Hp \ x ) 



(48) 



The boundary shift method discussed in this Section 
is a simple and straightforward approach to perturbation 
series expansions for dielectric cavities. As it is based on 
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Eqs. (|23[) - (f26|) . it first shrinks to zero the regions where 
the refractive index differs from its value for the circular 
cavity. Consequently, the calculation of the field distribu- 
tion in these regions remains unclear. Besides the direct 
continuation of perturbation series (|34[) inside these re- 
gions diverges. To clarify this point, we discuss in the 
next Section a different method without such drawback. 



Eq. (|49|) with 'potential' (|50|) can be rewritten in the form 

(A 2 + k 2 nl(x )) *(x) = -k 2 5n 2 (x )*(£ ) . (54) 

Then its formal solution is given by the following integral 
equation 



B. Green function method 

Fields in two-dimensional dielectric cavities obey the 
Helmholtz equations @ which can be written as one 
equation in the whole space for TM polarization 

(A + k 2 n 2 (x )) = (49) 

with position dependent 'potential' n 2 (x). For perturbed 
cavity (|20|l 

n 2 {x ) = nl(x ) +6n 2 (x) (50) 
where n^ix ) is the 'potential' for the pure circular cavity 

2r-\ _ j n2 when \x \ < R , . 

n ° [X) ~ { 1 when \x\ > R ' 

and the perturbation Sn 2 (x ) is equal to 

(n 2 - 1) when f(9) > and R < \x \ < R + \f(6) 
-{n 2 - 1) when f(0) < and R+ Xf(9) <\x\< R 
in all other cases (52) 

Hence, the integral of 5n 2 (x ) with an arbitrary function 
F(x) = F(r,6) can be calculated as follows 

, , r R+\f(8) 

/ Sn 2 (x )F(x )dx = (n 2 - 1) d6 F(r,6)rdr . 

(53) 



= -k 2 j G{x, y)6n 2 (y)V(y)dy (55) 

where G(x,y) is the Green function of the equation for 
the dielectric circular cavity which describes the field pro- 
duced at point x by the delta-function source situated at 
point y. The explicit expressions of this function are pre- 
sented in Appendix O 

It is convenient to divide the x plane into three circular 
regions r < Ri, R\ < r < i?2, and r > R 2 where 

Rl = mm e (R,R + Xf(9)) 

R 2 = maxe(R,R + \f(6)) . (56) 
The boundaries of these regions are indicated by dashed 
circles in Fig. |H Notice that the deformation 'potential' 
5n 2 (a? ) is nonzero only in the second region R% < r < R 2 . 
Due to singular character of the Green function (cf . Ap- 
pendix [C]), wave functions inside each region are repre- 
sented by different expressions. 

Let (r, 9) be the polar coordinates of point x. For sim- 
plicity we assume for a moment that f(9) < 0. Using 
(|C10jl . Eq. (|55p in the region r < R\ can be rewritten in 
the form 

*( 2 ) = Eto cos w4[*] (57) 



where £p[^] is the following integral operator 



x 2 (n 2 - 1) 



R 2 



>cos(p0) 



pdp 



_JjJpkri_ + i p _ H ii) {nx)J{knp)) 

_2Trxbp(x)Jp(nx) 4 1 v 



r 



(58) 



Assuming that we are looking for corrections to a quasi- 
stationary state of the non-perturbed circular cavity with 
the momentum equal m, one concludes that the quan- 
tized eigen-energies are fixed by the condition that the 
perturbation terms do not change zeroth order function 
(see e.g. [3]), i-e. 



t m m = i 



(59) 



which can be transformed into 

x 2 (n 2 -l) f 
S m (x) = ^ / d</>cos(p0) 



R+XfW 



*(p, <t>)pdp 



+ 



R 

iS m {x) 



J m {nkp) 



2irxJ m (nx) 
(H£> (knp)J m (nx) - H$(nx)J m (knp)) 



(60) 
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To perform the perturbation iteration of ([57)1 and ([60]) . 
integrals like the following must be calculated: 



V, 



1 



pm 



L p [J m (knp) cos(m</))] . (61) 



For small A the integral over p can be computed by ex- 
panding the integrand into a series of 8r = p — R 

R+^fW i 

F(p)dp » A/(^(i?) + -A 2 / 2 (0)F'(i?) + . • • • 

£ 1 

(62) 

Notice that this method is valid only outside the second 
region R\ < r < i? 2 which shrinks to zero when A — * 
(cf. (j56"|) ). In such manner, it leads to 



where 



and 



V pm = x 2 {n 2 - 1)(AV« + \ 2 V$) (63) 



v ® = ^M A ^ (64) 



V<$ = ^^[l + x(^^ + ^^)-2xS m (x)] 



* pm 



(65) 

Here A pm and B pm are defined in (|B5|) and (|B10[) . 
The second order terms can also be expressed through 

V, 



mp • 



J m (knr) 
J m (nx) 
■sr-^ J p (knr) 
Z-f J p (nx) 



cos(m8) 
cos(p8) 



(66) 



*pm 



The quantization condition (|59p in the second order 
states that 



i 



which can be expressed as 

S m (a;) = x(n 2 - l)(\A mm + ^A 2 [l + 



(67) 



(68) 



2x 



H^'(x) . 



Br, 



+ X 2 x 2 (n 2 -l) 2 



km 



k^rn 



Writing as in the previous Section x — xq + Xxi + X 2 X2 
where xq is a zero of S m (x) and using ()B1|) . one obtains 
the same series as ([33]) . Other expansions up to the sec- 
ond order also coincide with the ones presented in the 
Section HHSl 

To calculate the higher terms of the perturbation ex- 
pansion, the wave function must be known in the regions 
where the perturbation 'potential' 5n 2 (x) is non-zero. 
But exactly in these regions the Green function differs 



from the one used in Eq. (|58|) . In other words, a method 
must be found for the continuation of the wave functions 
defined in the first region r < Ri (or in the third one 
r > R2) into the second region i?i < r < i?2- 

The straightforward way of such a continuation is to 
use explicit formulae for the Green function in the sec- 
ond region and to perform the necessary calculations. 
As the radial derivative of the Green function is discon- 
tinuous, delta- function contributions will appear in cer- 
tain domains when calculating the integrals as in (|62|) . 
One can check that this singular contribution appears in 
the bulk only in the third order in A in agreement with 
Eqs. (P | and 

The expansion of wave functions into series of the 
Bessel functions (|57| . in general, diverges when r > Ri 
and the Green function method gives the correct continu- 
ation inside this region. Another equivalent method con- 
sists in a local expansion of wave function into power se- 
ries in small deviation from the boundary of convergence. 
As the value of the function and its radial derivative are 
assumed to be known along this boundary (r = Ri in 
our case) the knowledge of the wave equation inside and 
outside the cavity determines uniquely the wave function 
in the both regions. 



IV. APPLICABILITY OF PERTURBATION 
SERIES 

In the previous Section the formal construction of 
perturbation series has been performed for quasi-bound 
states in slightly deformed dielectric cavities. The pur- 
pose of this Section is to discuss in details the conditions 
of validity of such an expansion. 



From (|B5j) it follows that the coefficients A pm obey the 



inequality 



where £ stands for 



\A, 



e = A 



/ 


f(0) 




R 



< 2£ 



d6 : 



Sa 
^R2 



(69) 



(70) 



Here Sa is the surface where the perturbation 'potential' 
Sn 2 is non-zero and 7ri? 2 is the full area of the unper- 
turbed circle. The last equality is valid when (f2"Tj) is 
fulfilled which we always assume. 

Consequently, the perturbation formulae can be ap- 
plied providing 



£(n 2 - l)x Q ( 



1 



S p (x ) 



< 1 



(71) 



where (F p ) indicates the typical value of F p and x stands 
for Re(fci?) at a first approximation. 

The usual arguments to estimate this quantity for large 
Xo are the following. In the strict semiclassical approxi- 
mation, states with a corresponding incident angle larger 
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than the critical angle have a very small imaginary part 
and are practically true bound states. For closed circular 
cavities, the mean number of states (counting doublets 
only once) is given by the Weyl law 



N{Ej < n 2 k 2 ) 



An 2 
~8t7 



0{k) 



(72) 



where A = irR 2 is the full billiard area and n is the re- 
fractive index. The latter appears because by definition 



inside the cavity the energy is E 



,21,2 



For a dielec- 



tric circular cavity with radius R, the condition that the 
incidence angle is larger than the critical angle leads to 
the following effective area [l2| 



A eff(n) 



nR 2 i 



where 



- i R , 2 . R 

— - / (1 arcsm — )rdr 

Ir/u Ti" nr 

(•nx 

' 2 2 



2 

R 2 
4 



\fx 2 n 2 — m 2 dm 



(73) 

(74) 
(75) 



, 2 f ■ 1 1 L 1 
J ( arcsm — I — \j 1 — 



Here the first integral (|74|) corresponds to the straight- 
forward calculation of phase-space volume such that the 
incident angle is larger than the critical one and the sec- 
ond integral (Tf5|) is obtained from (fi"6|) taking into ac- 
count that x < m < nx. For n — 1.5, s n w 0.22. 

Consequently, the typical distance between two eigen- 
states is 



Sx 



(76) 



The eigen-momenta of non-confining eigenstates have 
imaginary parts of the order of unity (cf. Q13p ) and will 
be ignored. 

As S m (xo) = 0, we estimate that for p =/= m 



1 



1 



S'Sx 



Using (|Bip one finds that this value is of order of 



Op 

where constant C ~ 0.25n 2 s n /(n 2 — 1). 
leads to the conclusion that for typical Sj 
of applicability of perturbation series is 



(77) 



(78) 



s n — k 2 n 2 < 1 

87T 



With ||7I]> it 
the criterion 



(79) 



which up to a numerical factor agrees with 

But this statement is valid only in the mean. If there 
exist quasi-degeneracies of the non-perturbed spectrum 



(i.e. there exist p for which 1 / S p (x) is considerably larger 
that the estimate ([75])) then the standard perturbation 
treatment requires modifications. As circular cavities are 
integrable, the real parts of the strongly confined modes 
are statistically distributed as the Poisson sequences [l5[ 
and they do have a large number of quasi-degeneracies 
even for small k. For instance, these are double quasi- 
coincidences for the dielectric circular cavity with n = 1.5 



£14,2 


= 16.7170- 


0.03895 i, 


£11,4 


= 16.6976- 


0.4695 i 


£15,3 


= 17.5042- 


0.37540 i, 


£12,4 


= 17.5232- 


0.4612i 


^17,4 


= 21.5715- 


0.41621 i, 


£14,5 


= 21.5106- 


0.4712 i 


£25,1 


= 19.4799- 


0.00254 i, 


£21,2 


= 19.4830- 


0.1211 i 



We notice also a triple quasi-coincidence 

146,1 = 34.3110 - 2.2206 10" 6 i 
x 4h2 = 34.3167- 0.001982 i , 
£37,3 = 34.3 1 7 - 0.06408 i . 



(80) 



The existence of these quasi-degeneracies means that the 
perturbation series require modifications close to these 
values of kR for any small deformations of a circular cav- 
ity with n = 1.5 . 

Double quasi-degeneracy is the simplest case because 
there is only one eigenvalue of the dielectric circle, with 
quantum number, say p, which eigenvalue is close to the 
eigenvalue xq corresponding to the quantum number m. 
In such a situation, instead of the zeroth order equation 
(f3l"j) . the system of the following two equations must be 
considered 



S m (x (l + Sx))ai = M n ai 
S p (x (l + 5x))a 2 = M 2 \ai 



\- M 12 a 2 , 

\- M 22 a 2 (81) 

2 l)A ir 



where in the leading order My = xo{n" 

Expanding the S p functions and using the dominant 
order (|B 1 1) for S' p , the system (|5Tj) can be transformed 
into 



(si — 5x)ai 
(s 2 - 8x)a 2 



where 



si 



S m (xo) 

x (n 2 - 1) 



-A n 



S2 



-4i2a2 , 

^21^1 

xo(n 2 



1) 



- A 



22 



(82) 



(83) 



Our usual choice is 5 m (xo) = but for symmetry we do 
not impose it. The compatibility of the system ([82ll leads 
to the equation 



Sx — Si Ayi 

A 2 \ 5x - s 2 



= 



(84) 



Its solution which tends to s± when A\ 2 A 2 \ — > is 



A l t U l ( \ A i 4A 12^21 

ox = - si + s 2 + -(si - S2)4/l + 7 

2 2 V [S1-S2V 



= Si 



A 12 A 21 



•S'2 



si 1 + a/1 + AA 12 A 21 /( Sl - s 2 ) 2 



(85) 



10 



e 




FIG. 5: Cut disk cavity. eR is the distance between the cut 
and the circular boundary. 



When A12A21/ {s\ — S2) 2 is small, Sx in (|85|) tends to the 
usual contribution of the second order ([33]) . Therefore 
in this approximation, expression (|33[) may be used for 
all p except the one which is quasi-degenerate with Xq. 
For this later one, the whole expression ([85| (without 
Si) has to be used. A useful approximation proposed 
in [lj| consists in taking the modification (f8"5)) for all 
the non-degenerate levels which reduces the numerical 
calculations. 

As the circular billiard is integrable, the probability 
of having three and more quasi-degeneracies is not neg- 
ligible (cf. ([80|0 . The necessary modifications can be 
performed for any number of levels but the resulting for- 
mulae become cumbersome. In Appendix |D] we present 
the formulae for three quasi-degenerate levels. 



V. CUT DISK 

As a specific example, we consider a deformation of 
the circular cavity which is useful for experimental and 
technological points of view [171] . Namely a circle is cut 
over a straight line (see Fig. [5]). Such a deformation is 
characterized by the parameter e <C 1 which determines 
the distance from the cut to the circular boundary. This 
shape corresponds to the following choice of the defor- 
mation function f(9) 



56» 4 
~24 



(86) 



when 16*1 < 9 m and f{6) =0 for other values of 9. Here 
6 m is the small angle 



6 m = arccos(l — e) 



7 2e- 



V2 
12 ' 



3/2 



(87) 



Using the formulae discussed in the preceding Sections, 
we compute all the necessary quantities and compare 
them with the results of the direct numerical simulations 
based on a boundary element representation similar to 
the one discussed in 11311 - 
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FIG. 6: Comparison between perturbation and numerical 
spectra for the cut disk with e = 0.05. Open black circles 
are results of the direct numerical simulations and red full 
circles correspond to the perturbation expansion (I33I) . 



The spectrum of quasi-bound states for a cut disk with 
e = 0.05 is plotted in Fig. [6] To get a good agreement 
in the region Ke(kR) ~ 15 — 20, it is necessary to take 
into account double quasi-degeneracies and, in the region 
close to Re(fe-R) = 35, triple degeneracy ([50)) has been 
considered. The agreement is quite good even in the 
region of large ki? where many families intersect. 

Two wave functions of this cut disk are plotted in 
Fig. [7] The first one is obtained by direct numerical sim- 
ulations and the second one corresponds to the pertur- 
bation expansion. Even tiny details are well reproduced 
by perturbation computations. In the direct numerical 
simulations, wave functions are reconstructed from the 
knowledge of the boundary currents. This procedure re- 
quires the integration of the Hankel function, Hq 1 ^ [l3[. 
As this function has a logarithmic singularity, a small re- 
gion around the cavity boundary has been removed to 
reduce numerical errors. It explains the white region 
in Fig. [7^). For an easier comparison, the same region 
has been removed also from the perturbation result in 
Fig. [7b). 

Finally, in Fig. [SJ the far-field emission pattern com- 
puted numerically for the cut disk is compared with the 
same deduced from the perturbation series (|48|) . The 
both are normalized to unit maximum. Once more a 
good agreement is found. The approximate positions of 
the main peaks in the far-field pattern correspond to the 
diffracted rays emanated from the discontinuities of the 
cut disk boundaries and reflected at the critical angle ([5]) 
on the circular boundary (see Fig. [5]). If (0%, 82, ■ ■ ■) and 
(9'i, #2) • • •) are the directions of two such refracted rays 
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* s 



< » 



a) b) 

FIG. 7: a) Modulus square of the wave function for the cut disk with e = .05 obtained by direct numerical simulations 
corresponding to quasi-stationary eigen-momentum kR = 16.655 — 0.0199 i b) The same but calculated within perturbation 
expansion with m = 21 and p = 0. The corresponding eigen-momentum is kR = 16.659 — 0.0191 i. The high intensity regions 
are indicated in black. 
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FIG. 8: Far field emission pattern, |F((9)| 2 , for the same state 
as in Fig. [7] Results from direct numerical simulations are 
plotted in solid black line, while the dashed red line indicates 
the perturbation result (I48|) . The blue points correspond to 
experimental results with e = 0.05 and R — 60 /im. 

one gets from geometrical considerations 

9 2 = A9 C - 6 m - | , 9' 2 = A9 C + 9 m - | . (88) 

Here 9 C is the critical angle © and 9 m is defined in (|8"T|) . 
For e = .05 and n = 1.5, 9 Y w 155.4°, 9[ w 168.2°, 




FIG. 9: Dashed and dotted lines: two main diffracted rays 
responsible for dominant peaks in the far-field pattern of a 
dielectric cut disk. All rays hit the circular boundary wi th 
an angle equal the critical angle l[8|. 



6» 2 w 59°, and 9' 2 w 95.4° which agrees with Fig-U 

To complete this study, far-field experiments have been 
carried out according to the set-up described in [l8| . The 
cavities are made of a layer of polymethylmethacrylate 
(PMMA) doped by 4 — dicyanomethylene — 2 — methyl — 
6 — (4 — dimethylaminostyryl) — AH — pyran (DCM, 5 
% in weight) on a silica on silicon wafer. To obtain a 
good resolution of the shape even for such a small cut as 
e = 0.05, cavities are defined with electron beam lithog- 
raphy (Leica EBPG 5000+) by C. Ulysse (Laboratoire 
de Photonique et de Nanostructures, CNRS-UPR20). A 
scanning electron microscope image of such a cavity is 
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FIG. 10: Scanning electron microscope image of a cut disk 
etched with an electron beam. 

shown in Fig. |TD] As specified in [l8j], the cavities are 
uniformly pumped one by one from the top at 532 nm 
with a pulsed doubled Nd:YAG laser. The light emit- 
ted from the cavity is collected in its plane with a lens 
leading to a 10° apex angle aperture. The directions of 
emission are symmetrical about the 0° axis according to 
the obvious symmetry of the cut-disk shape. In Fig. [5J 
the intensity detected in the far-field is plotted versus the 
angle with blue points. The position of the maximal 
peak around 160° is reproducible from cavity to cavity 
with e = 0.05 and agrees with both numerical and per- 
turbation approaches. 

The example of the cut disk clearly demonstrates the 
usefulness of the perturbation method presented in this 
paper for deformed circular cavities. 

VI. CONCLUSION 

We considered in details the construction of perturba- 
tion series for deformed dielectric circular cavities. The 
obtained formulae can be applied for the calculation of 
the spectrum and the wave functions as well as other 
characteristics of these dielectric cavities (e.g. far-field 
emission patterns). We checked these formulae on the 
example of the cut disk which is of interest from an ex- 
perimental point of view. Cavities of other shapes (e.g. 
spiral [H) 

can be considered analogously. This method 
can also be used to calculate the influence of a small 
boundary roughness on the emission properties of circu- 
lar cavities. 
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APPENDIX A: WHISPERING GALLERY MODES 

The purpose of this Appendix is to calculate the 
asymptotic of quasi-stationary eigenvalues for a fixed ra- 
dial quantum number p and a large azimuthal quantum 
number m — > oo. 

These resonances are well confined, so let us first 
find the corresponding asymptotic expression for zeros 
of Bessel functions J m (x m ) = 0. To achieve this task it 
is convenient to use Langer's formulae (see e.g. Q 7.13.4) 
which are valid with 0(m~ 4 / 3 ) accuracy 

J m (x) = T(w)[J 1/3 (z)cos(7r/6) - Yi/ 3 (z)sin(7r/6)] , 
Y m (x) = TH[J 1/3 (z)cos(7r/6) + y 1/3 (z)sin(7r/6)] , 

where T(w) — w^ 1 / 2 (w — arctan^)) 1 / 2 , w — (x 2 /m 2 — 
l) 1 / 2 , and z = m(w — arctan(w)) . 

The combination of the Bessel functions is expressed 
as follows 

Ji/3{z) cos(tt/6) - Y 1/3 (z) sin(7r/6) 

= 3 1 /6 2 V3 z -i/3 Ai ( _ ( 3 z/2) 2/3 ) (Al) 

where Ai(x) is the Airy function 

1 f°° t 3 

Ai(x) = - / cos (— + tx)dt . (A2) 
it Jo 3 

Therefore in the intermediate region when z is fixed and 
m 3> 1, the zeros of the J m Bessel function correspond 
to 

z = \)T ( A3 ) 

where rj p are the modulus of the zeros of the Airy func- 
tion Ai(— 77 P ) = (rji re 2.338 , tj 2 ~ 4.088 , % ps 5.521). 
Finally, the Bessel function zeros have the following ex- 
pansion 

x m = m + am 1 / 3 + /3m~ 1/3 + C(to _1 ) (A4) 

where 

a = 2- 1 '\, ^-L^y 2 . (A5) 

For p = 1 this expression agrees numerically with the one 
given in 0] sect. 7.9. 

The whispering gallery zeros of Bessel functions (|A4|) 
permit to calculate explicitly the whispering gallery 
modes for the dielectric disk. Indeed we are interesting 
in the solutions of equation 

l-f(nx) = (A6) 

in the region xn close to m, what means 

xn = x m + Sx (A7) 
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with Sx <C x m . The expansion of the right-hand side of 
Eq. (fM} leads to 



H {1) ' x 

11 m / ^ 



H m n 



i 2 ' 3 y/n 2 — 1 



(A8) 

In the left-hand side of Eq. (|A6[) , both numerator and 
denominator can be expanded into powers of Sx taking 
into account that J m {x m ) = 

/' /' 4- St J" 4- (St) 2 V" 

J m {Xm + dX) ^ SxJ^ + (Sx) 2 J^/2 + (Sx)^/6 {Xm) ■ 

(A9) 

Using the Bessel equation 

JM + - z J' m {z) + (l - J m (z) = (A10) 

one can check that all but one derivatives Jm^ / J' m {x rn ) 
are at most ©(to -1 ) and can be neglected. The only 
exception is 

J 1 " 2a 

fW = 273 ) • ( A11 ) 

J' m m z/A 
Finally expansion (|A9|) takes the form 

^±( Xm + Sx) = ^--Sx-^+Oim- 1 ) . (A12) 
J m ox 3to/' j 

Combining this equation with (|A6|) and (|A8|) one obtains 

3 2 



5 a; 



n an ^ 

v ^T^l + m 2/Z v ( n 2 Z 1)3/2 (p _1 ) + °( m ) - 

(A13) 

These formulae lead to Eq. (JT8| . 

Other 'whispering gallery' modes correspond to x close 
to to. With the same notations as above, Langer's for- 
mula can be written for the Hankel function 

H£>(x) = T{w)e™' & H { ^{z) + 0(m- 4 / 3 ) . (A14) 

3 (z) exis 
7.11. 



The zeros of H 1 / 3 (z) exist only when z = re 17r with real 
r. From formula Ml 7.11.42, it follows that 



SUITTU 

+ c -^^^Hi 2 Hr) . (A15) 
sm7w 



In this way we get 



x [J 1/3 (r) cos(7r/6)-y 1/3 (r)sin(7r/6)] . (A16) 

This is the same combination of the Bessel functions as 
Eq. IAip for J m [x). Therefore the first complex zeros of 

the Hankel function, Hm\xm)i have a form similar to 

H 

x m = m + &m 1/3 + /3to~ 1/3 + C(to -1 ) (A17) 



where a = e _27 "/ 3 a and /3 = e _4,n / 3 /3 with the same a 
and /? as in (|A5p . 

The next step is to find the asymptotic of J m (nx) for 
complex x = x rn . As 

Im(f m ) = -2 -1/3 sin(27r/3)?7pTO 1/3 -|-0(TO -1/3 ) , (A18) 

it tends to — oo with increasing of to. Therefore from 
(Tj"4")) it follows that, instead of the cos(. . .) term, only the 
positive exponent has to be taken into account. Then the 
left-hand side of (IA6D leads to 



V 



(nx m ) = u 1- — — + G(to- 1 ) (A19) 



i^Hll 



(n 2 - 1)to 2 / 3 



+ C(to- 1 ) . 



As H$(z) al so ob eys Eq. (jATOj) one obtains the same 
expansion as (|A12() 

H \r(x m + St)) = -I - Sx^jt + (Dim- 1 ) . (A20) 



Sx 3to 2 / 3 



Eqs. (|A19|) and (|A20[) lead to the following value of Sx 
with C(to _1 ) accuracy 

l-^V (A21) 



Sx = — - 



\/?i 2 — 1 m 2 / 3 \A?-T V 3 
Collecting these equations we get (fT9l . 



APPENDIX B: CONSTRUCTION OF 
PERTURBATION SERIES 

From Eqs. (f2"4")) - the following relations are valid 
at the circle r — R: 

*i-*a = -A 2 ^ 7p cos(^) + 0(A 3 ) 
p 

9*2 



<9r <9r 



fe [5m (a;) cos(to#) 
A (a p + Xf3p)S p (x) cos(pO) 

p^m 

R (t)> 



a 2 *i 9 2 * 2 

<9r 2 <9r 2 



1? 



[(5 m (x) + x(n 2 - 1)) cos(m6>) 



A ^(^(x) + x(n 2 - 1)) cos(p6>) 

p^m 



o(\ 2 



and 



9 3 *i 9 3 * 2 



Q r 3 Q r 3 



kSm (x) 



m 2 + 2 



R 2 



; 2 2 

k n 



H-M' „2i 

-fc 3 (r^- l)^U-(x) +fc 2 

' J m 



i? 



cos(m<2) + 0(A) . 
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In the zeroth order S m (xo) = 0. Expanding S m (x) with 
x as in ([50)) into a series in A, one gets 



dS 

S m (x + Xxi + X 2 x 2 ) = Xxi-^- 

-<D{\ 2 ) 



> / dS m 1 2® 2 S, 



2 1 da; 2 



where all derivatives of S m are taken at x = xq. These 
derivatives are deduced from the Bessel equation (|A10I) 

^L(x)=-(n 2 -l)--S m (x) 
ox X 

(1)' \ 



S m {x) ( S m (x) + 2 *" (x) 

\ Hm i 



In particular, when x = xq 



dx 



■(x Q ) = -K-i) 



(Bl) 



d 2 s., 



1 



(1)' 



In the A first order, it leads to 

— x\(t? — 1) cos(m0) + a p 5 p (a;o) cos(p6>) 

/(*) 



-a;o(n — 1) cos(m6 



Coefficients a p and the first eigenvalue correction, x±, 
are determined by comparison of the Fourier harmonics 
in both sides of Eq. (|B2l 



and 



(B3) 



(B4) 



where A pm are the Fourier harmonics of the deformation 
function 

A pm = ^- ( f(6) cos(p0) cos(m0)d0 . (B5) 

TT-fl JO 



Here 



2 for p f 
1 for p = 



(B6) 



In the A second order, it leads to the following two equa- 
tions 

]T 7p cos(p#) = \x 2 {n 2 - l)^cos(m0) (B7) 



(B2) and 



/ 1 ( H [1)l 



R 



dS, 



H 



(i)' 



+ 7>7 



H, 



(i)' 



H, 



(i) 



008(7716*) 



cos(p6>) 



a p (^S p + x (n 2 - lfj cos(p9) 



p^m 



(1)1- 



— -x {n - 1) I 1 - x — ^- cos( ? 



(B8) 



Using (|B2|) the right hand side of equation (|B8|) can be 
rewritten as 



. From (|B7|) it follows that for all p 



(n 2 — 1) I icq ^ cos(p0) + x\ cos(m9) 



p^m 



r(l)' N 



f (9) I Hri 



lp = 2 X o( n ~ l ) B pm ■ 



(B9) 



where B pm represents the Fourier harmonics of the 
square of the deformation function 



Unknown coefficients can be determined by equating the 
Fourier harmonics in both parts of equations (|B7p and 



f 2 (9)cos( P 9)cos(m9)d9 . (BIO) 
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For p 7^ m Eq. (|B8|) gives 

as, 



h, 



i)' 



f3 p S p + xiotp-^- + 7 P — ^yy = (n 2 - 1) xiAp, 



r(l)' 



+x ^ a fe A pfc + 2 X o(^l +x '' ( ^ jBp m 

The mfi 1 harmonic of the same equation determines the 
second correction to the quasi-stationary eigenvalue 



1 2A ,o H ™'\ 



(1)' 



(n 2 - 1)H, 



(i) 



-£o(l + 2^0 

2 v ff, 



(1)' % 

I B mm 



xq u k A mk - x\A r , 

k^rn 



Rearranging these equations and using the first order re- 
sults and Eq. (|B9j) . one gets 



X 2 = Xq 

n, 



(Bll) 



and 



/3p = x - 



S n 



-^ B P m I 1 + X °{ 







pm A mrn 


/Xq dSp 

\S P dx 


rr(l)/ 


"'")) 


o-(l) + 


4" ') 



(B12) 



k^r. 



APPENDIX C: GREEN FUNCTION FOR THE 
DIELECTRIC CIRCULAR CAVITY 

The Green function of the dielectric Helmholtz equa- 
tion for the circular cavity, G(x,y), is defined as the so- 
lution of the following equation 

{A s + nt(x)k 2 )G(x,y) = 6(x-y) (CI) 

where is the 'potential for the circular cavity defined 

in (EE]). 

Let us first consider the case when the y source point is 
inside the circle. In this case, when the x point is inside 
the cavity, the advanced Green function has the form 

00 

G(x,y) = V A m J m (nkr)e im ^ 



m= — oo 



H^(kn\x-y\ 



and when the x point is outside the circle 

00 

G(x,y) = B m Hg\kr)e" n ^ 



(C2) 



(C3) 



Here and below we assume that points x and y have polar 
coordinates (r, 9) and (p, 4>) respectively. 

Constants A m and B m are calculated from the bound- 
ary conditions on the interface using the expansion Q 

7.15.29 of H^\k\x -y\) 

OO 

H^(k\x-y\)= J m (kr)H^\k P y m ^ (C4) 

m— — oo 

when r < p and 

00 

H^(k\x-y\) = H^(kr)J m (kp)e im ^ (C5) 



m— — 00 



when r > p. 

For the TM polarization the Green function and its 
normal derivative at circle boundary are continuous. 
That leads to the following system of equations (with 
x = kR) 

A m J m (nx) + — J m {nkp)H${nx) = B m H$(x) 
4i 

nA m J' m (nx) + J J m (nkp)H^'(nx) = B m H^'(x). 
Its solutions are 

A m ■ 

B m ■■ 



H£\x) . H£\nx) 

J m [knp) - — -J m [knp), 



2i:xA m J m (nx) A\J m (nx) 
J m (knp) . 



2ttxA t 



where 



A m = nJ' m (nx)Hg)(x)- J m (nx)H$'(x) (C6) 
= J m {nx)H^> (x)S m (x) . 

In deriving these expressions, the Wronskian (0 7.11.29) 



has been used 



2i 



J u (x)HW(x) - J' v {x)H™{x) = — . (C7) 

TTX 

For the y source point outside the circle, when x is inside 
the cavity, then 

00 

G(x,y)= J2 C m J m (nkr)e im ^ (C8) 

in— — oo 

and when x is outside the circle, then 

00 

G(x,y) = D m H^(kr)e im ^ (C9) 



m=— oo 

+ L H U(k\x-y\) . 
Constants C m and D m are computed exactly as A m and 



Dm = J (nx) g (i )(M _^(£L g (i) (M 



4iH m (x)~ 
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The final expressions of the Green function follow from 
the above formulae. 

When the y source point is inside the circle, the plane 
is divided into three parts: 1) r < p, 2) p < r < R, 3) 
R < r. Denoting the Green function in these regions by 
the corresponding numbers, it can be written as 



APPENDIX D: THREE DEGENERATE LEVELS 



G j (x,y) = ^gW( r ,p)co S \p(6 -<(>)} 



p=0 



(CIO) 



where 



(1) J p (nkr)J p (nkp) 
9 P {r,p)~ 2 ^e p j2( nx j Sp ^ 

4i Jp(nx) ' * 

(2) / \ _ _J_ J P (nkr)J p (nkp) 
9 P KT,P)~ 2jrx £ P j2( nx)Sp ( x) 

4^TTT^ 1) ( hf ) J pH " H^(nx)J p (knr)], 

t:1 tj rp\TtX J 



f\r,p) = - 



1 J p (nkp)H^\kr) 
2^ P J p {nx)H { p 1] {x)S p {x) 



where e p was defined in (|B6j) . 

When point y is outside the circle, the plane is divided 
into three different regions 1) r < R, 2) R < r < p, 
3) p < r. With the same notation as above, the Green 
function can be written 



Gj (x, y) = J2 ~9p } ^ P) C0S lP( e ~ ft] ( cn ) 

where 



For three quasi-degenerate levels, instead of Eq. ([84] 
one gets the 3x3 determinant 



Sx-sx A 12 A 13 
A n Sx-s 2 A 23 =0 (Dl) 
A 31 A 32 Sx - s 3 



which leads to the cubic equation 

(fa) 3 - o-i (fa) 2 + (o- 2 - a)Sx - a 3 + (3 = . (D2) 
Here <jj are the elementary symmetric functions of Si 
crx = s 1 +s 2 + s 3 , a 2 = sis 2 + s 2 s 3 + s 3 si , a 3 = S1S2S3 ■ 
and (because Aij = Aji) 



a = A\ 



A 2 

2 ^ ^23 



A 2 
^31 ) 



(3 = 2A 12 A 23 A 31 + Sl A 2 23 + s 2 A 2 31 + s 3 A\ 2 . 



To solve Eq. (|D2|1 . the next steps are standard. The 
substitution 



1 

fa = V + gCTl 

transforms Eq. (|D2[) to the reduced form 



- PV + q = 



(D3) 



(D4) 



where 

P = - 



1 2 1 

-a 2 +a 2 -a , q = -— o\ + -ai<7 2 ~<7 3 + 5 . (D5) 



with 5 = (3 — \<Jia. Finally after the transformation 



y 



3z 



(D6) 



\ 1 J p (nkr)H p 1] (kp) 



2^ J p (nx)H i p 1) (x)S p (x) ' 

9P [ ,P) ~ 2itx H^{x)S p {x) 
+ T^^m^K ] ( x )Mkr) - HW(kr)J p (x)} , 



4i " 



~9 { p\riP) = --^-e p 



1 Hf\krW(kp) 



2nx Hp(x)S p (x) 

+ ^P^ L J^[Hp [1 \x)Jp{kp) - H£\kp)J p {x)] . 
41 H p '(x) 

Notice that Gx(x_,y) = G 2 (y,x), G 3 (x,y) = G x {y,x), 
and G 2 (x,y) — G 3 (y,x). It means that in all cases the 
Green function is symmetric: G(x,y) — G(y,x) as it 
should be (see e.g. [1J|). 



one gets the equation z 3 — ^fpr + q = which is a 
quadratic equation in variable 



Its solution is 



10 



w 



± 



p 
27 



(D7) 



(D8) 



Eqs. (|D3[) . (|D6|) . (|D7|) . and (|D8J) give the well known 
solution of the cubic equation (|D2[) . The question is how 
to choose a branch which tends to si when — > 0. The 
discriminant of this equation is 



4 



27 



cf- 



(D9) 



where 



1 



27 



-G\a 2 - a 3 



If 

27 V 



ct 2 



-en 
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and 



X r2 l z( 2 3 1 

3 i 1 2/ ^ 2\ ^ / 1 2\2 

- 7fj a +g" (0-2- 3 CT l)- gQt(CT2- gO-l) 



Using the identity 



(si - s 2 )(s 2 - s 3 )(s 3 - Sl) 



the expression (|D8j) can be transformed into 

w = w - -S + 7^( s i ~ S 2)(S2 - S 3 )(s 3 - si) 



'1- 



108e 



[Ol - S 2 )(s 2 ~ S 3 )(s 3 - Sl)] 2 



Finally the root of the cubic equation (|D2p which tends 
to si when Aij — > is 



z = —s 
3 



, 27(5 3^/31. 

2^ + ^ (Sl " S2)(S2 " S3)(S3 " Sl) 



1/3 



108e 



[(si - s 2 )(s 2 - s 3 )(s 3 - si)] 2 



where 



1 3 1 1 
wo = ~ 01 ~ g°"i02 + 2 CT 3 



+ — ^=(S1 - S 2 )(S 2 - S 3 )(S 3 -Si) 

= ^(s 1+ s 2 e- 2 -/ 3 + S3 e 2 -/ 3 ) 3 



where s = si + s 2 e 27ri / 3 + s 3 e 27r1 / 3 . 
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